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Abstract 

The hierarchy of fermion masses in the standard model may arise via the breaking of discrete 
gauge symmetries. The renormalizable interactions of the flavor-symmetry-breaking potential can 
have accidental global symmetries that are spontaneously broken, leading to pseudo-goldstone 
bosons that may drive inflation. We consider two-field, axion-monodromy inflation models in 
which the inflaton is identified with a linear combination of pseudo-goldstone bosons of the flavor 
sector. We show that the resulting models are nontrivially constrained by current cosmological 
data as well as the requirements of viable flavor model building. 
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I. INTRODUCTION 


The prevailing approach to solving the horizon and flatness problems of conventional 
Big Bang cosmology is inflation, a period in which the nniverse underwent exponential 
expansion due to the effects of the nearly constant energy density provided by a scalar 
field [1]. Models of inflation are often studied in terms of the properties of the inflaton 
potential, with somewhat less focus on other roles the inflaton might play in extensions of 
the standard model. If the inflaton has no purpose other than to provide the source of the 
energy density that drives inflation, then model building becomes isomorphic to studying 
ways of generating different functional forms for the inflaton potential. These possibilities, 
now cataloged (see for example H), differ in their detailed predictions for the spectrum of 
fluctuations in the microwave background that are observed in experiments like Planck j3] 
and BICEP2 |1]. 

In this paper, we consider a scenario in which the inflaton is an integral component of 
an extension of the standard model that aims to address one of its substantial mysteries: 
the hierarchy of elementary fermion masses. Models of flavor based on horizontal discrete 
symmetries postulate that these symmetries are broken via a set of fields, called flavons, 
that couple to standard model fermions through higher-dimension operators. Discrete flavor 
symmetries can often lead to accidental continuous global symmetries among the renormal- 
izable terms of the flavon potential. In the present work, we consider the possibility that the 
inflaton may be identified as a linear combination of the approximate goldstone bosons that 
arise when these accidental symmetries are spontaneously broken. We will ultimately be 
interested in two-field models of inflation, for reasons described below, which distinguishes 
the present work from the relatively sparse literature that explores the use of flavon fields 
for a similar purpose [5]. 

Consider the simplest possibility, a Z^v flavor symmetry under which a single flavon field 
$ transforms as <I> —)■ Ci;<h, where u = exp(27ri/iV). If the fermions of the standard model are 
charged under the discrete group, then a tree-level Yukawa coupling that would otherwise be 
forbidden can arise via a higher-dimension operator. For example, for a down-type quark, 
one might have 

^Q^HcjyWn + h.c. , ( 1 . 1 ) 

where H is the standard model Higgs doublet, Mp is the flavor scale, and p is an integer. The 


2 



Yukawa coupling is associated with the ratio {{(j))/MFy which can be much less than one; 
operators with different values of p can easily provide a hierarchical pattern of entries in the 
associated Yukawa matrix. If less than the Planck scale, the scale Mp is identihed with that 
of new heavy states that account for the origin of the higher-dimension operators. However, 
a simpler assumption, that we adopt henceforth, is that the scale Mp is the reduced Planck 
scale M*; the desired operators appear as part of the most general set that are allowed 
by the local symmetries of the theory, as one expects based on our current understanding 
of quantum gravity [B]. An immediate implication of our assumptions is that the vacuum 
expectation value (vev) ($) < M*, which will provide an important constraint in our attempt 
to identify the inflaton with a part of the held <h. 

To obtain an inhaton potential that is sufficiently hat, we require that the goldstone 
boson degree-of-freedom receives no contributions to its potential from renormalizable terms 
involving $. Let us therefore assume that N > 5. The renormalizable terms in the potential 
are simply 

I/($) =+. (1.2) 


Terms such as (<I>^ -|- h.c.) are forbidden by the Z^v symmetry. Using the nonlinear decom¬ 
position 


$ = 


^ + f 
V2 


exp{te/f) , 


(1.3) 


where f/V2 = ($), one sees immediately that U(<I>) is independent of 6, he., the potential 
V{9) is exactly hat. The potential in Eq. ( 1.2[ ) has an accidental global U(l) symmetry 
and the held 6 is the goldstone boson that results from its spontaneous breaking. Global 
symmetries are not respected by quantum gravitational corrections, so it is no surprise that 
there are Planck suppressed corrections. 


that generate a potential for 6, where cq is an unknown order-one coefficient. Planck- 
suppressed operators that directly break the discrete havor symmetry are not present since 
we assume in this example (and will require in all our models henceforth) that we work with 
discrete gauge symmetries, which satisfy appropriate anomaly cancellation conditions and 
are immune to quantum gravitational corrections. For the reader who is unfamiliar with 
discrete gauge symmetries, we review the basic issues relevant to our model building in the 
appendix. 
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The operator in Eq. (1.4) leads to the 9 potential 


rw = c„ M." [1 - cos (/ve//)| , (1.5) 

where we have added a constant so that V (0) = 0. This is nothing more than the potential of 
“Natnral Inflation” scenarios |7]. However, this potential is not adequate for our purposes. 
It is well known that if one requires that Natural Inflation provides ~ 50 — 60 e-folds of 
inflation and predicts a spectral index within the range allowed by current Planck data, 
then / must be well above the Planck scale [8] . For our present application, this would imply 
that (0)/M* is not a small flavor-symmetry-breaking parameter and we lose the ability to 
predict standard model Yukawa couplings in a controlled approximation. 

We therefore must consider other ways of generating potentials for the pseudo-goldstone 
inflaton that allow a sub-Planckian decay content /. The options assuming a single held 
inhation model are limited. For example, models of “multi-natural” inhation [S], in which one 
arranges for additional sinusoidal terms in the potential, can accommodate a sub-Planckian 
havon vev, but tend to predict = 0.95 in this limit [9], at the very edge of the 95% 
exclusion region following from Planck data. A different class of models that can more 
easily provide cosmological predictions consistent with Planck data are two-held models of 
the axion monodromy type [TUHTBj . We will show that these can be adapted for the present 
purpose. 

The two pseudo-goldstone helds can have their origin if there are two havon helds, $ and 
X, that transform under the discrete group Z* x Z^. We assume that each held transforms 
only under one of the Zjy factors, 

$ and x ^xX , (1-6) 

where = exp(27rz/p) and = exp(27rz/r), where p and r are integers. For p > 5 and 
r > 5, the renormalizable terms in the potential are 

I/($, x) = - m^x^x + ^(X^X)' + A^^t^x^x , (1.7) 

where Xp is a portal-type coupling. There is no difficulty in choosing parameters such that 
each held develops a vev. This potential has an accidental U(l)xU(l) global symmetry that 
is spontaneously broken. Extending our previous parameterization, we write 

® exp{i9/fg) and x = exp(ip//p) . (1.8) 
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Spontaneous symmetry breaking renders the fields 0o and xo massive so that they are 
decoupled from the inflation dynamics. The potential for the goldstone bosons V{p, 9) that 


follows from Eq. (1.7) is exactly fiat. 

We will discuss later how to generate a potential for p and 9 of the following axion- 
monodromy form 


V{p,9)=k\ 


1 + cos ( Y 

Jp 


+ A 


np 9 

1 - cos I —-— 

Jp Je 


(1.9) 


where n is an integer. The first few terms in the expansion of the first cosine factor have 
the same form as —m^r^/2 + Ar.r^/4!, the shift-symmetry-breaking potential W{r) assumed 
in the Dante’s Waterfall scenario discussed in Ref. na. In that work, W{r) was assumed 
to be generated by non-perturbative effects associated with moduli stabilization in string 
theory, as for example in Ref. m- In this paper, we only consider field theoretic origins 


of the potential, where the emergence of the functional form given in Eq. (1.9) is readily 
obtained. For the purposes of graphical display, if one plots the potential as if p and 9 were 
polar coordinates, one would find a “hill” generated by the first cosine factor, circumscribed 
by a descending spiral “trench” generated by the second. Inflationary trajectories track the 
minimum of the trench. As 9 advances by 271 fp along the trench, the p coordinate does not 
return to the same value; this monodromy allows for large numbers of e-folds to be achieved 
within a bounded, sub-Planckian region of field space. We assume that the decay constant 
fg satisfies 


—= = XM^ Tn 0.22 , 

V2 


( 1 . 10 ) 


where A is a fiavor-symmetry-breaking parameter of the same size as the Cabibbo angle. 
This will allow us to identify the field $ (and perhaps in some models both $ and x) 
fiavons that can be used in flavor model building. We will see that the discrete symmetry 
Zp X Zjf serves four purposes: (i) it assures that there are goldstone bosons that have no 
potential generated by renormalizable couplings, {ii) it will serve as a flavor symmetry to 
create a hierarchy of standard model fermion Yukawa couplings, {Hi) it will lead to the 
correct pattern of couplings in a new gauge sector that provides for the desired form of the 


infiaton potential, Eq. (1.9), and (iv) it will keep quantum gravitational corrections to the 
potential highly suppressed. 

Our paper is organized as follows. In the next section, we discuss the inflationary dy¬ 


namics that follows from the potential given in Eq. (1.9). We identify solutions in which 
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inflation ends when single-field slow-roll conditions are violated and other solutions where 


the termination of inflation is analogous to a hybrid model HU. In Sec. |IIH we consider 
model building issues, in particular, how the discrete symmetries of the theory play an im¬ 
portant role in assuring that we obtain the proper potential, and how the same symmetries 
can be used to produce a plausible model of standard model fermion masses. In the hnal 
section, we summarize our conclusions. A brief appendix is provided to review relevant facts 
about discrete gauge symmetries. 


II. INFLATONARY TRAJECTORIES 


In this section, we consider inflationary trajectories in the two-field potential given by 


Eq. (1.9), that are compatible with flavor model-building requirement Eq. (1.10). We give 


two example solutions that differ qualitatively in how inflation ends. 


Termination without a waterfall. 


For our Erst solution, we make the parameter choice fp = fg = fi and also define 


fi/n = f 2 - We assume /i 3> / 2 , which is equivalent to n S> 1. The potential Eq. (1.9) then 
takes the form 


V{p,d)=K\ 


1 cos I Y 


+ K 


1 \ p ^ 

1 — cos --^ 


( 2 . 1 ) 


/2 h 

The second cosine term creates a series of trenches on the surface of the potential defined 
by the first cosine term. If the held 9 is plotted as a polar coordinate, the trenches form 


spirals originating at p = 0. As in Ref. 
basis p = cp + s6 and 9 = c9 — sp with 

/i 


, it is convenient to work in the rotated held 


c = 


This allows us to rewrite the potential as 


and 


s = 


h 




( 2 . 2 ) 


V(p,ff) = At 


1 -|- COS 


cp -|- s 9 


+ A 


. P 

COSI- 


(2.3) 
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where / = /i/ 2 /\/ fi + /I- The modulations in the potential due to the cos(p//) term 
create the trench, whose location is given by dV/dp = 0, or 

cp + s6 


p\ At . 

sm 1 -J - seisin 


h 


= 0 


(2.4) 


The inflaton is the linear combination of the helds that slowly rolls along the trench; inflation 
terminates when the slow-roll conditions are violated. For the solutions considered in this 
subsection, the stability condition d’^V/dp^ > 0 will hold throughout this trajectory. 

To study inflationary observables, we first consider a good approximation to the single- 
held inhaton potential, which holds for our choice of parameters and can be studied analyt¬ 
ically, and then discuss an exact numerical approach that we use to conhrm the validity of 
our results. Let us dehne n = sc{h.\/h. 2 ) and consider parameter choices where k 1. It 


follows from Eq. (2.4) that to good approximation 


p// S3 2irj , 


(2.5) 


where j is an integer. Given our assumption that /i ^ / 2 , it follows from Eqs. (2.3)-(2.5) 


that d‘^V{p, 6)/dp^ > 0, conhrming that the trench is stable. Substituting Eq. (2.5) into our 
original potential yields 

V{e)= At\l +cos (6 + 9/fo)] , ( 2 . 6 ) 


where 6 = 2ti sej and /o = fi/s. Setting j = 0 is equivalent to redehning the origin of 
held space, so we will ignore 6 henceforth. We note that the present approximation scheme 
dihers from the one used in Ref. [12], in which one would expand the sinusoidal functions 


in Eq. (2.4) to linear order in their arguments, but is nonetheless accurate as we conhrm 


numerically later. We note that s 1 in the limit u S> 1, so that the derived quantity /o 
can be super-Planckian even when the decay constants fi and /2 are not. 

We compare the predictions of the model to the latest results from the Planck Collabo¬ 
ration [3]. The slow roll parameters are dehned by 


1 

^ “ 2 


R' 


R" , R'R'" 

rj = — and 7 = 


(2.7) 


V / ’ ' V ' 1/2 

where the primes refer to derivatives with respect to the inhaton held and we work in units 
where the reduced Planck mass M* = Mp/\/+n = 1. In the present case, these are given by 

1 

VI 




( 2 . 8 ) 
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V = 


cos{9/fo) 


(2.9) 

( 2 . 10 ) 


/o 1 + cos(0//o) 

7 =-4tan2[^/(2/o)] • 

Jo 

Inflation ends when e{6f) = 1. The initial value of the inflaton, 6i is determined by the 
requirement that we achieve a desired number of e-folds of inflation, given in general by 


N = 


X 




d9 = 2 /q In 


sin[gy/(2/o)] 

sin[0i/(2/o)] 


( 2 . 11 ) 


We set = 60 in the numerical results that follow. We evaluate the slow-roll parameters 
and the potential V{9) at 9i in determining the spectral index n* = 1 — 6 e -|- 2r], the ratio of 
tensor-to-scalar amplitudes r = 16e, the running of the spectral index = 16er] — 24e^ — 27 


and the scalar amplitude Aj^ = l//(247r^e). From Eqs. (2.8)-(2.10), it follows that 


ns = l + ^(l-2sec%/{2fo)]) , (2.12) 

Jo ^ ^ 

r =-^taia^[9i/{2fo)] , (2.13) 

Jo 

=-^tan^[0i/(2/o)]sec^[0i/(2/o)] , (2.14) 

Jo 

(l + cos[0i//o]) csc'^[6i/fo] • (2.15) 

To illustrate a viable solution, consider the parameter choice (again, in units where M* = 1) 

/i = 0 . 22^2 , (2.16) 

/2 = / 1 / 2 I , (2.17) 

= A 2 = 0.006 , (2.18) 


which corresponds to n = 21 and k ~ 1/21. We find that the initial and final fields for the 
inflaton trajectory are given by 

(p, 9)i = (6.04 X 10-^ 6.74) and (p, 9)f = (1.50 x 10-^ 19.14) , (2.19) 

respectively. Using this value for 9i, we hnd the following set of cosmological parameters: 


Us = 0.96, (2.20) 

r = 0.060, ( 2 . 21 ) 

= -0.00046, (2.22) 

= 2.2 X 10“® . (2.23) 
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p(x 10^) 


FIG. 1: Path followed by the inflaton during 60 e-folds of inflation corresponding to the 
Sec. IIA in units where M* = 1. The background is a density plot where darker zones 
values of the potential than lighter ones. 


solution of 
have lower 


Fig. displays the path followed by the inflaton during the 60 e-folds of inflation for this 


particular solution. The predictions in Eq. (2.23) are consistent with the results from the 


Planck experiment [3]: = 0.968 ± 0.006, r < 0.12 (95% C.L.), Ur = —0.003 ± 0.007 and 

A|, = 2.19 ± 0.08 X 10“®. (The value of A|j, also called Ag, was taken from the hrst column 
of Table 3 in Ref. [3].) 

We may check the validity of the results in this section by numerically evaluating the 
slow-roll parameters in the two-held problem. Let a represent the linear combination of the 


helds that evolves along the minimum of the trench. Given that da = y dp^ + dO'^ along the 
trench, it follows that we can write the derivative of the potential with respect to a as 


dJ^V 

da^ 


^ ^ dp\ d 


dO 


tr 


dO 


v[erp{e)tr 


(2.24) 


where the subscript “tr” indicates quantities evaluated along p{0)tr, the solution to Eq. (2.4) 


Note that as the quantity da is dehned above, the kinetic terms for a are canonically nor¬ 


malized. The slow roll parameters can be evaluated numerically according to Eq. (2.24). 
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We find in this case that Us = 0.96, r = 0.060, Ur = —0.00046 and = 2.2 x 10 in 


agreement with the results in Eq. (2.23). 


B. Termination with a waterfall. 


For different choices of the model parameters, inflation will end before e = 1 is reached, at 


a point where there is no longer a solution to Eq. (2.4). At this point, the stability condition 
d'^V/dp^ > 0 is also not satisfied, and the fields evolve rapidly in a direction orthogonal 
to the original trajectory If one visualizes the motion by plotting the fields as polar 
coordinates, the evolution corresponds to a transition from spiraling to rapid motion in the 
radial direction, eventually ending at a global minimum. In Ref. [12] this was called the 
waterfall, in analogy to the behavior of hybrid inflation models ini, where stability in one 
field direction can be a function of the value of a second field. 

Given an input of model parameters, we determine the final inflaton field value by 
solving 

r)2T/ 

= 0 , (2.25) 


dp^ 


tr 


and then the initial value a,- from 


N = 


r°'f 

V 

J ai 

V' 


da 


(2.26) 


where the primes refer to derivatives evaluated numerically according to Eq. (2.24), and a 


(~ 6) is the canonically normalized inflaton field. Again, we set N = 60. To illustrate a 
solution that ends with the waterfall behavior, consider the parameter choices 


/i = 0.22^2 , 

/2 = /i/17 , 

Ai = 3.38 X 10“^ 
As = 1.61 X 10“^ 


(2.27) 

(2.28) 

(2.29) 

(2.30) 


which corresponds to n = 17 and k, = 1.13. We find that the initial and final fields for the 
inflaton trajectory are given by 


(p, e)i = (6.83 X 10■^ 1.63) and (p, 0)/ = (0.0281, 5.2970) , 


(2.31) 
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respectively. Using this value for 6i, we find the following set of cosmological parameters: 


Us = 0.96, 
r = 0.0078, 

Ur = -7.2 X 10"^ 


= 2.2 X 10"^ 


(2.32) 

(2.33) 

(2.34) 

(2.35) 


These are consistent with the ranges allowed by Planck, as quoted in the previous subsection. 


The complete inflaton trajectory, extending beyond the point where Eq. (2.4) is no longer 
satished, can be found by solving the coupled equations of motion 


dV 

p + 3Hp+-^ = 0 , 
. dV 

^ + = 0 > 


(2.36) 


where H is the Hubble parameter. The result is shown in Fig. One can see from the plot 


that the bottom of the trench given by Eq. (2.4), denoted by the thick red line, approximates 


the actual trajectory, given by the thin green line, very well. The inflaton eventually oscillates 
about and then settles at the global minimum of the potential. 

III. MODELS 

A. Origin of the potential 


The successful inflation potentials presented in the previous section correspond to a po¬ 
tential of the form 

e 


v{p,e) = A 


1 + cos ( ^ 
Ip 


+ A 


cos 


np 

T 


(3.1) 


where n is an integer. Here we consider the possibility that this potential arises via the 
effects of anomalies associated with new gauge groups. 

Hence, we extend the standard model gauge group by the additional factors 


SU(iVi)xSU(A 2 ), and introduce the fermions Al ~ Ar ~ (Ni, 1) and B^l’ ~ B'p’ Cl ^ 


tf ~ B® 

Cr (1,N2). We would like the Lagrangian to contain the following interactions 


C D H^ArAlX + hfB^^Bfx + h^CRCL^* + h 


c. 


(3.2) 


2=1 
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FIG. 2: Inflaton trajectory, in p-0 space, overlaid on a contour plot of the potential, in units where 
M* = 1. The bottom of the trench is indicated by the thick red line while the inflation trajectory 
is denoted by the thin green line. 


Here, the are Yukawa couplings and the terms shown generate heavy fermion masses 
when the $ and y fields develop vevs. The accidental global U(l) symmetries are each 
chiral when appropriate charges are assigned to the H, B and C fermions. However, these 
symmetries are anomalous with respect to the new gauge groups. Triangle diagrams lead to 
the interactions [13] 


9l 


327r2 V/, 


^ 1 FiFi + 


9l 

327r2 


'f-UF.F. 

Ip JO, 


(3.3) 


and the non-perturbative generation of a potential 


V{p,d)=A\ 


1 I P 
1 - cos I — 

Jp 


+ A 


np e 

1 - cos I --— 

Jp JO 


(3.4) 


with the scales Ai and A 2 identihed with the scale of strong dynamics for each SU{N) factor. 
(We assume Ni and W are chosen so that each group is asymptotically free.) Redehning 
the origin of held space via 


p ^ p + vr/p and 6^6 + nnfe 


(3.5) 


puts the potential in the form that we previously assumed in Eq. (3.1). Note that the new 
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gauge groups may be spontaneously broken at a scale well below Ai and A 2 without affecting 
our conclusions. 


The interactions given in Eq. (3.2) are clearly not generic. In the absence of our discrete 


charge assignments for $ and y, there would be no reason for the $ field to avoid coupling to 
the A and 5-type fermions directly, nor would there be any prohibition of explicit fermion 
mass terms. Hence, this sector is suggestive of additional symmetries even had we not 
put them forward immediately as a starting assumption in our model building. Given the 


transformation properties of $ and y fields under the x symmetry, Eq. (1.6), we can 


account for the desired pattern on couplings in Eq. (|3.2|) by choosing 

Ar —)■ UJx^R y 




Cl —)■ u^iCl , 


(3.6) 


with the remaining heavy fermions taken as singlets under the discrete group. However, we 
must now enlarge the fermion content to assure that discrete gauge anomalies are cancelled 
(see the appendix), and do so in a way that assures that the additional fermions can become 
massive. To demonstrate that this can be accomplished, let us consider an example suggested 


by one of our previous cosmological solutions, discussed in Sec. IIA corresponding to the 


potential in Eq. (3.1) with n = 21. Let us choose p = r = 21. First, we note that there are 
21 5-type fermions transforming each with Z 21 charge -|-1, where we specify the charge Q by 
defining the group element to be exp(2i7rQ/21). This implies that the Z 2 i-SU(A^" 2 )^ discrete 
anomaly cancellation condition would be satisfied by the 5 particle content alone. The A 
and C fermions, on the other hand, lead to anomalies, so we include additional fermions 
with matching gauge quantum numbers and the discrete transformation rules 


A® ^ 


,10 /id) 


(t = 1...2) 


Lxl Ul 


, 10 ^ 6 ) 


cf -»eg’ 


('i = 1...2) 


(3.7) 


which allow the anomaly cancellation conditions to be satisfied. Finally, we note that 
these fields will develop masses as a result of Planck-suppressed operators, for exam¬ 
ple, -I- h.c. and /M^ + h.c. , which lead to masses of order 

~ 10^1 GeV. 


The discrete symmetry that we have assumed to assure the form of couplings in Eq. (3.2) 


also leads to a suppression of direct Planck suppressed corrections to the potential. Since 
quantum gravitational effects must respect the discrete gauge symmetry, the lowest order 


13 








Qil 

Q 2 L 

Qsl 



^R 


^R 

b'k 

6 

5 

3 

2 

-1 

-3 

-1 

-2 

-2 


Lil 

L 2 L 

Lsl 



^R 

^IR 

^2R 

^3R 

0 

0 

0 

5 

3 

1 

-3 

-3 

-3 


TABLE I: charge assignments q, where the group transformation is defined by exp(2i7rg/21). 

The Higgs doublet is a singlet under the flavor symmetry. 

operators that will correct the potential have the form or the scale of 

these corrections are of order ~ negligible compared to the values of Ai 

and A 2 that we found previously to be of order 

B. Standard Model Flavor 

The helds $ and x can now be utilized in constructing models of standard model fermion 
masses. These helds will appear in higher-dimension operators that generate the small entries 
of the standard model Yukawa matrices. Given our choice = (x)/M* = A, the size 

of these entries will be determined by powers of the Cabibbo angle A. In this subsection, we 
present one example in which the desired set of higher-dimension operators is obtained via 
the same discrete symmetries that were used to obtain the inhaton potential. We focus on 
the n = p = r = 21 model just discussed, in which the <I> and x fields each transform under 
a separate Z 21 symmetry. Of course, other choices of the symmetry group are possible, and 
the present choice does not suggest a unique set of fermion charge assignments (since there 
are many possible Yukawa textures that are viable). The example we give here will suffice 
by serving as a proof of principle^. 

The simplest incorporation of the n = 21 model in a flavor sector is via the identification 
of Zfj as the flavor symmetry and $ as the sole fiavon field. The charge assignments of the 
standard model fermions and a set of right-handed neutrinos are given in Table Entries 

^ It should also be clear that one could alternatively construct a model starting with the n = 17 potential 
that we discussed earlier, but there are no new qualitative insights gained by presenting two very similar 
examples. 
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of the Yukawa matrices arise from Zf^-invariant higher dimension operators. For example, 
the 1-1 entry in the up-sector Yukawa matrix involves the helds Qh^Hur, which has flavor 
charge —8. This arises at lowest order via 

+ h.c. , (3.8) 

and hence the corresponding Yukawa matrix entry is of order A®. Since cj® and are 
identical, there is another possible operator, + h.c., but it is of higher 

order and can be neglected. We may populate the remaining entries of the quark and charged 
lepton Yukawa matrices in a similar manner. We hnd 



^A® A^ A®^ 


^ A^ A^^ 


^ A® A® A ^ 


w = 

A^ A" A2 

, Yd = 

A" A® A® 

, n = 

A® A® A 

, (3.9) 


^A® A^ l) 


^ A2 A X j 


^ A® A® A ^ 



where order one coefficients in each entry have been suppressed. These achieve the desired 
ratios rriu/mt ~ A®, rric/mt ~ A"^, md/rrii, ~ A^, and ms/rrib ~ A^, with the charged lepton 
Yukawa mass eigenvalues comparable in size to those of the down quark sector. It is not 
hard to verify that the choice of right-handed neutrino charge assignments leads via the 
see-saw mechanism to a neutrino mass matrix of the form /Ar\Y^, where Ar is the 
right-handed neutrino mass scale, (H) is the standard model Higgs vev, and is a matrix in 
which each entry is of order A° times a function of (typically many) undetermined order one 
coefficients. These can be chosen to obtain the desired phenomenology without unnaturally 
large or small values of the individual coefficients^. 

Finally, we must check that the standard model fermion charge assignments in this model 
satisfy the linear Ibahez-Ross anomaly cancellation conditions for the non-Abelian gauge 
groups and gravity. Summing the Zf) charges times the appropriate multiplicity factors for 
the color SU(3), weak SU(2), and gravitational anomalies gives 21, 42 and 63, respectively. 
These results mod 21 are zero, indicating that the discrete gauge anomaly cancellation 
conditions discussed in Appendix remain satished. 

^ It is not necessarily the case that an alternative model that predicts the neutrino mass hierarchy via 
powers of A is more desirable than this example. The reason is that the predictions for neutrino mass 
matrix entries in such a model also come multiplied by functions of products of a number of the order 
one operator coefficients. This can spoil the naive A power counting without any individual operator 
coefficient being unnaturally small or large. This is a problem that is unique to the neutrino sector in 
such models when the mass matrix arises via the seesaw mechanism. 
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IV. CONCLUSIONS 


Models of standard model flavor that are based on discrete gange symmetries can have 
accidental continnons global symmetries that are spontaneonsly broken. We have argned 
that a linear combination of the approximate goldstone bosons that may arise in these 
models can serve plansibly as the inflaton in two-held models of inhation based on the 
axion monodromy idea. These models can accommodate the cnrrent Planck data on the 
microwave backgronnd [3] while allowing the havor-symmetry-breaking vacnnm expectation 
valnes (vevs) to remain snb-Planckian. This is important in the present work since the ra¬ 
tios of the havon vevs to the rednced Planck scale serve as small havor-symmetry-breaking 
parameters in onr models, which allows one to predict the standard model Yukawa coupling 
entries in a controlled approximation. In addition to making correct Yukawa coupling pre¬ 
dictions possible, the discrete symmetries of the theory also maintain the correct pattern of 
the interactions in a new gauge sector, leading to the desired form of the inhaton potential; 
they also keep the quantum gravitational corrections to the potential well under control. 
The literature on models of standard model fermion masses is vast and it is imaginable that 
more economical and compelling examples of havor-sector inhation models are yet to be 
found. The present work suggests that exploring the full landscape of such models may be 
a fertile direction for future investigation. 
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Appendix A: Discrete Gauge Symmetries, Briefly 

It is well known that continuous gauge symmetries are not violated by quantum gravita¬ 
tional ehects. Under what circumstances is the same true for discrete symmetries? It was 
noted long ago by Ibanez and Ross (IR) [H] that a discrete group that arises as a subgroup 
of a continuous gauge symmetry inherits this protection. While the full theory must satisfy 
the anomaly cancellation conditions relevant for the continuous gauge groups, IR determined 
the conditions that are relevant in the low-energy theory, below the scale at which the con- 
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tinuous gauge symmetries are broken. Since some of the fermions in the complete theory 
may become massive and decoupled when symmetry breaking occurs, the low-energy theory 
includes only part of the fermion content that contributes to anomaly cancellation in the 
full theory. The low-energy constraints should refer only to the light fermion content, which 
in the present context corresponds to models dehned below the reduced Planck scale M*. 
If the appropriate consistency conditions are satished, the discrete gauge symmetry can be 
treated as fundamental, without reference to specihc high-energy embeddings. 

The constraints that we apply in our model building are the linear IR conditions in¬ 
volving non-Abelian gauge group factors; these follow from triangle diagrams involving two 
non-Abelian gauge group factors and one factor of the continuous gauge group in which 
the discrete symmetry is embedded. For example, the ZAr-SU(M)^ anomaly cancellation 
condition is [18] 

l^C.9. = lrAf. (Al) 


Here r is an integer, g* is the Z^v charge of the fermion (which transforms under Zat by 
exp[i27rqi/N]) and Ci is the Casimir invariant given by Tr(T“T^) = where the T“ are 

SU(M) generators in the representation of the fermion. Since all the fermions in the 


model presented in Sec. Ill are in the fundamental representations of the relevant SU(M) 
gauge groups, Ci = 1/2; the linear IR conditions simply requires that the Zat charges of 
the fermions that transform under a specihed SU(M) factor sum to an integer multiple of 


N. According to IR, when N is odd (relevant to the model of Sec. Ill) the gravitational 
anomalies linear in Zat are cancelled when the sum of all the Zat charges are also an integer 
multiple of N. It is straightforward to verify that these conditions are satished by the charge 
assignments displayed in Table [1} 

What about the other possible anomaly cancellation conditions? First, IR note that the 
linear conditions involving the Abelian gauge groups do not lead to any useful constraints on 
the low-energy theory [T8|. Banks and Dine (BD) [T9| later showed that the IR conditions 
non-linear in the discrete group make a tacit assumption about the high-energy embedding 
of the theory, through the requirement that both the light and the heavy fermions have 
integer U(l) charges. BD show that there are consistent, non-anomalous theories (ones in 
which the effective discrete symmetry at low energies is smaller than that of the full theory) 
in which the low-energy spectrum does not satisfy the non-linear IR constraints; their failure 


17 




only implies the existence of heavy fermions with fractional charges. Thus, the non-linear 
IR conditions are not required for the consistency of the low-energy effective theory. BD 
note that the surviving discrete anomaly cancellation conditions are physically sensible: for 
example, the condition for the cancellation of the ZAr-SU(M)^ anomaly also guarantees that 
there are no t’Hooft interactions generated by SU(M) instantons that would explicitly break 
the ZTV symmetry. This physical constraint pni is completely independent of the high-energy 
embedding. 
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